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EFFECT OF THE LIFT COEFFICIENT ON PROPELLER FLUTTER 
By Theodore Theodoreen &nd Arthur A. Regier 

SUMMARY 


Flutter of propellers at high angles of attack Is 
discussed, and flutter data obtained In connection with 
tests of models of large wind-tunnel propellers are 
analyzed and results presented. It Is shown that In the 
high angle-of-attack range flutter of a propeller lnvarla 
bly occurs at a speed substantially below the classical 
flutter speed. The angle of attack at which flutter 
occurs appears to be nearly constant and independent of 
the initial blade setting. Thus, the blade simply twists 
to the critical position and flutter starts. Formulas 
have been developed which give an operating angle in 
terms of the design angle and other associated parameters 
and these relations are presented in the form of graphs. 
It is seen that the flutter speed is lowered as the 
initial design lift coefficient is increased. It Is 
further shown that by use of a proper camber of the pro- 
peller section the flutter speed may approach the classi- 
cal value. A camber for which the blade will not twist 
Is found to exist, and the corresponding lift coefficient 
Is shown to be of special significance. 
r s Xzit 

The bias si cal )f lutter speed and the divergence speed 
of a propeller are shown to be approximately the same 
because of the centrifugal-force effects. It appears 
that a propeller will not flutter until the blade twists 
to a stall condition near the divergence speed. Strobo- 
scopic observation of several propellers confirmed this 
theory. It was observed that, regardless of the initial 
pitch setting of the propeller, the blades always twisted 
to the stall condition before flutter commenced. The 
problem of predicting propeller flutter is thus resolved 
primarily into the calculation of the speed at which the 
propeller will stall. 
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INTRODUCTION 


The present study of propeller flutter was conducted 
In connection with the design of several large wind-tunnel 
propellers for the Langley, Ames, and Cleveland Labora- 
tories of the NACA. Vuind- tunnel propellers are not 
required to operate In a fully stalled condition and can 
therefore be designed with small margin of safety against 
flutter. Airplane propellers, on the other hand, must 
have a considerable margin of safety since they are 
required to operate in the stall or near-stall condition 
in take-off. The results of the present tests are of 
wide interest since they apply to the general problem of 
the effect of high lift coefficients on the flutter 
velocity of a propeller. • 

There are two principal types of flutters (1) "clas- 
sical” flutter and (2) "stall” flutter. Classical flutter 
is an oscillatory instability of an airfoil operating in 
a potential flow. The problem of classical flutter was 
solved theoretically in reference 1. Stall flutter 
Involves separation of the flow and occurs on airfoils 
operating near or in the stall condition of flow. Studer 
(reference 2) studied this type of flutter experimentally 
with an airfoil in two-dimensional flow. He found that 
the stall flutter speed was very much lower than the 
classical flutter speed and that, as the angle of attack 
of the airfoil was increased, the change from classical 
flutter to stall flutter was rather abrupt. 

The problem of propeller flutter Is somewhat dif- 
ferent from the problem of wing flutter In that the 
change between classical flutter and stall flutter appears 
to be much more gradual. This gradual change of flutter 
speed with angle of attack has not been clearly under- 
stood, and attempts to calculate the flutter speed of 
propellers operating under normal loads have not been 
entirely satisfactory. 

A third type of flutter, which may be referred to as 
”wake” flutter because it occurs on propellers operating 
at zero lift In their own wake, has sometimes been 
observed. Self- excited torsional oscillations of the 
propeller blade occur at frequencies which are integral 
multiples of the rotational speed of the propeller. At 
low speed the frequency. of oscillation is equal to the 
torsional frequency of the propeller blade In still air. 
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The first such oscillation appears when the propeller 
rotational speed reaches approximately one-twentieth of 
the blade torsional frequency. This oscillation dis- 
appears as the rotational speed is increased but reap- 
pears at each integral multiple of the propeller speed 
until the classical flutter speed is reached. This type 
of flutter is not important for normal propeller opera- 
tion. 

SYMBOLS 

L representative length of propeller blade 

c chord of propeller section 

b chord 

t tl.ioKr.ess of propeller section 

R radius to tip of propeller 

r radius to propeller section 

K torsional stiffness of representative section 

q dvns.'r.ic pressure of relative air stream 

P density 

k ratio of mass of cylinder of air of diameter 

equal to chord of airfoil to mass of airfoil 

a angle of attack 

A a angle of twist or deformation of blade at repre 

sentative section 

o mo angle of attack for which there is no twist 

a l Q angle of attack for zero lift 

C mc ^ moment coefficient about quarter- chord point 

Cl lift coefficient 
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°Lu 


a ^ 


w a 

“h 

M 

V D 

v f 

v f. 


untwisted or design value of Cl 

lift coefficient for ideal no-twist condition ■ 

location of center of '"gravity as measured from 
leading edge 

location of center of gravity with reference to 
elastic axis in terms o'f semi chord 

coordinate of torsional stiffness axis in terms 
of semichord as measured from mid chord posi- 
tion , 

nondimenslonal. radius of gyration of airfoil 
section in terms of semi chord referred to a 

torsional frequency, radians per second 
> • 

bending frequency, radians' per second 

Mach number 

divergence speed 

flutter speed 

flutter speed corrected for compressibility 


Subscripts : 


u 

cr 

I 

c 


untwisted or design 
critical 


i 


ideal 

compressible 

incompressible 
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EFFECT OF LOADING ON THE TWIST OF A PROPELLER BLADE 


The centrifugal force on a propeller has a component 
in the direction perpendicular to the relative flow, which 
is very nearly equal to the aerodynamic force. This 
statement is exactly true if the propeller is designed to 
avoid Pending stresses and is approximately true in any 
case, since the bending forces are small compared with 
the aerodynamic forces. (See fig. 1.) The twist of the 
propeller at some representative section may be expressed 
by the relation 

Aa K = qLc 2 -^x - -0-^ (<*u + Aa - (1) 

The value of the critical velocity q^ for which 

divergence occurs is obtained from equation (1). Diver- 
gence evidently occurs for the condition 

Act oo 


or 


Aa 

au + Aa - a mo 




For Aa—^oo then. 



( 2 ) 


By substitution of this value of q cr for q In equa- 
tion (1) the twist becomes 


Aa 



%>) 


L cr 


' 1 - 


L cr 


( 3 ) 
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The moment coefficient around the quarter- chord point 
may he written as 

°“cA = ( ai o ‘ <tao)( x ' 

This equation may he rewritten as 



The relation between the untwisted or design value of 
the lift coefficient C]^ and the actual or measured 

value of Cl resulting from the twist may he expressed 
as 

dC L 

= -aST (“» - ai o + 


dC L 

Ct — Ct + Ag 

L Hi da 

By substitution for Aa from equation (3) and by use of 

the value of !% from equation ( i|_) the following equa- 

o 

tion is obtained for Cl: 


Aa 


) 


(5) 



q. 

^cr 


t 1 -. 


<L_ 

cr 


x - 
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or the fallowing equivalent relation is obtained: 


C Lu = °L - ^ (°L + 


( 6 ) 


The increase In lift coefficient due to twist Is evidently 

°*B/k 


40L = ^( 0L+ 1 


X - 


(7) 


There is no increase in Cl, or no twist, if 


'La 


-- 

[ 1 
X “ 4 


( 8 ) 


where C^ indicates the value of Cl at which no 

twist occurs. This relation Is plotted In figure 2. 

This figure shows that, for the Clark Y airfoil with 
center of gravity at 1 j 4 percent and = -0.07, the 


value of the lift coefficient at which no twist of the 
blade is incurred Is 0.37* In this case the angle of 
zero twist Is not very far from the Ideal angle of attack 
of the Clark Y airfoil, which is about O.i^.0. The following 
discussion shows that It is desirable to operate the pro- 
peller at the Ideal angle of attack since operation at 
this angle delays stall and thus obviously causes an 
Increase in the flutter speed. (See reference 3 for 
discussion of ideal angle of attack.) 


A propeller. If generated as a true helix, will not 
be subjected to any centrifugal twisting moment; In fact, 
if the blade width of the propeller is adjusted to achieve 
the desired blade loading at an angle everywhere propor- 
tional to the helix angle, there will be no twist. This 
statement must be modified slightly, however, because the 
radial generating lines through the leading and trailing 
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edges may not be continued to the center and, as a result, 
the angle at the- tip may be decreased. By proper plan 
form and distribution of mass, therefore, a propeller 
may be designed to have zero twist. A highly tapered 
propeller tends to decrease its tip angle due to cen- 
trifugal twisting moment, and a more nearly rectangular 
plan form induces an increase in the tip angle. These 
effects are in reality small compared with the indirect 
effect due to the aerodynamic forces resulting from the 
bending of the blade. 


EXPERIMENTAL STUDIES OF FLUTTER OF PROPELLERS 
AT HIGH LOADINGS 


A number of propellers of different designs, some 
representing existing NACa wind-tunnel propellers for 
which data were available and others representing pro- 
posed wind-tunnel propellers, were tested as wind-tunnel 
fans in a small open tunnel. A cross-sectional sketch 
of the test setup is shown In figure 3 . 

The lift coefficient of the blades was changed by 
changing the area of the tunnel exit. The value of C]^ 

was calculated from the relative wind direction at the 
0.8-radius station and .the angle of attack of the 
untwisted blade. 

The propeller tips were observed by stroboscope 
through a small window in the tunnel wall, centered In 
the plane of the propeller. Blade-tip deflection and 
twist and the pressure Increase of the air passing 
through the fan were recorded. These data were used to 
give an independent check on the operating lift coeffi- 
cient. Owing to the nonuniformity of the blades, the 
variation of the lift coefficient with the radius, and 
other causes, the value of the operating lift coefficient' 
could be determined only within about 0.1 to 0.2. 

Although a number of propellers were tested, only 
results from the tests of two fairly representative pro- 
pellers are reported. The propellers were made of lami- 
nated spruce and had flat-bottom Clark Y sections. Pro- 
peller A, for which data are given In figure ij., was a 
six-blade propeller 45 inches In diameter. Propeller B 
was a single-blade propeller having the same diameter but 
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two-thirds the chord and thickness of propeller A. The 
propeller of smaller. oro.ss saotlon, propeller . B, -was used 
to reduce the flutter speed and to make possible a study 
of the flutter modes. Propeller B was tested in the same 
location in the tunnel as the other propellers, but a 
booster fan was attached at the rear of the motor to 
force the air through the tunnel during the tests of pro- 
peller B. 

The vibration frequencies of the propellers are as 
follows: 


Mode 

Vibration frequency 
(cps) 

Propeller A 

Propeller B 

First bending 

1 

714- 

I 4 I 4 . 

| Second bending 

21+6 

172 

Torsion 

555 

314-0 


Torsion and bending strain gages were attached to 
propeller B, and the flutter amplitudes and frequencies 
were recorded. The results of these tests are given in 
figure 5 . The flutter speed was changed by changing the 
blade lift coefficient. It may be observed that the 
bending amplitude is large and the flutter frequency is 
low ( 17 O ops) at the highest observed flutter speed. In 
a lower range of flutter speed - from l+OO to 500 feet per 
second - the observed bending amplitude is small and the 
flutter frequency attains a higher value, 2 lj.O cycles per 
second Instead of 170. This flutter evidently involves 
the second bending mode, whereas the flutter at top speed 
Involves the first bending mode. At the lowest flutter 
speed and the highest blade loading, the flutter reverts 
to a condition of pure torsional oscillation at a fre- 
quency corresponding to that of pure torsion measured In 
free air, namely, 3^0 cycles per second. 

The results of the flutter tests of propeller A are 
shown in figure 6 , which Is a plot of equation ( 6 ) for 
Ci^ = 0.37* the value of Cl^ for a Clark Y flat- 

bottom airfoil of 12-percent thickness. The straight 

lines that converge at -3— = 1.0 and Cl,. = 0.37 in, 

<*cr 
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figure 6 are lines representing constant values of Cl* 
Figure 6 shows that, if the design lift coefficient Ci^ 
is 0 . 6 , for instance, the blade will twist so that Cl 

is 1.0 at -S— = O. 63 . Diagrams similar to figure 6 may 
^cr 

be used to obtain Cl for any propeller in terms of Ci^. 

Figure 6 also shows curves for data plotted with 

uncorreoted for compressibility and for the same data 
plotted with CLu corrected for compressibility by 
Glauert's formula 


C Lu 


c 


C lTi! 

Vl - 


The data corrected for compressibility show that, when 
the propeller was set with the blade at stall, the flutter 
speed q/q cr was only O.17. As was decreased 

from O.85 to O.65, q/q cr increased from G.37 to 0.62. 

It may be noted that in this range the flutter occurred 
at an approximately constant Cl of 1.1. A further 
increase in q/q cr caused a rather sharp drop in Cj, 

for flutter. The data uncorrected for compressibility 
show that, as Ci^ varies from. O.7S to 0 . 55 » Gl ^ or 

flutter varies from 1.1 to 0.82. The validity of the 
compressibility correction as applied to Cl in figure 6 
may be questioned, but results of tho tests clearly 
indicate that the propeller twisted to a lift coefficient 
near unity before flutter occurred. The propeller can 
carry a lift coefficient exceeding unity without flutter 
If q/q cr is low enough; however, as the classical 

flutter speed is approached, the flutter lift coefficient 
becomes less than unity. In other words, the amount of 
stall necessary to excite flutter is less as the classical 
flutter speed Is approached. This conclusion Is In agree- 
ment with the experiment of reference 2 and seems logical 
in view of the various flutter modes described by figure 5 

The minimum flutter 3 peed for a completely stalled 
propeller Is of importance for propellers that operate 
in such a condition at times - in take-off, for example. 
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There is no reliable method known for calculating this 
minimum flutter speed for a completely, sjtalled propeller, 
and further Investigation of the problem is required. ' As 
the angle of attack is increased from a normal value, the 
flutter speed will decrease until a minimum is reached. 
There / is evidence that this m inimum is related to the 
von Karman vortex street; the minimum probably corresponds 
to a coincidence of the torsional frequency with that of 
the von Karm&n vortex street, in. figure- 6 - it may be 
observed that this minimum appears at . — c 0.17 for 

^cr 

propeller A. This value may be fair for wooden propellers 
but cannot be taken as valid for metal propellers. Flut- 
ter on metal propellers operating In the stall condition 
has been observed to occur at values of q/q cr as low 

as 0.0i|.. 

Figure 7 is of interest as a verification of the 
theoretical treatment in this paper. It gives the twist 
of the propeller tip as a function of Ci^ for a constant 

propeller speed at = O.37. This twist was observed 

‘Icr 

by means of telescope and stroboscope. The line on the 
figure is drawn through the point for an angle of twist 
of 0° and a value of Ci^ of 0.37* a3 predicted by fig- 
ure 2 . The slope of this line, which may be calculated 
by use of equation (3)» is adjusted to fit the data. 
Equation ( 3 ) is based on a simplified propeller and gives 
the twist at the representative section. At Ci^ = O.78 

and = 0 . 37 » equation (3) gives Aa = 2 .ij.°. The 

^cr 

observed twist at the propeller tip for this condition 
was 5.1° . It was observed by direct measurement that the 
torsional stiffness at the representative section waa 
2.3 times that at the tip. The observed tip twist is 
therefore consistent with the expected v&lue. 


DETERMINATION OF DIVERGENCE SPEED AND 
CLASSICAL FLUTTER SPEED 


The. divergence speed of a propeller can be calculated 
from equation ( 2 ) if proper values are selected for L, 
c, and K but may be more conveniently found from the 
formula of reference ij. (p. 17) which Is 
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(9) 


Reference Ij. (p. 17) also gives an approximate flutter 
formula that appears to hold very well for a heavy wing 
and for small values of Oh/wa " the conditions for 
normal propellers. The approximate flutter formula is 


v f 



( 10 ) 


where x a is the location of the center of gravity with 
reference to the elastic axis in terras of the semichord. 

It may he noted that formulas (9) and (10) are alike 
except for an additional term x a in formula (10). If 
the center-of -gravity location coincides with the location 
of the elastic axis, x a = 0, ..the two formulas are iden- 
tical. 


It was shown In the discussion of figure 1 that 
there are two moments acting about the elastic axis, the 
aerodynamic-force moment and the centrifugal-force moment. 
If the aerodynamic forces are balanced by components of 
the centrifugal force, the resulting moment Is the same 
as If the aerodynamic forces acted with a moment taken 
about the center of gravity of the airfoil section. For 
prooon *.-.re , therefore, the dynamic- stiffness axis may be 
tai-.or. at c .incur of gravity and the divergence speed 

of the prop}licr will be given approximately by the flut- 
ter formula, equation (10), or vj) = vf. It may be men- 
tioned that for normal propellers the location of the 
elastic axis and the location of the center of gravity 
are usually very close together. 

The location of the center of gravity frcm equa- 
tion (10), a + x a , is expressed in terms of the semi- 
chord as measured from the midchord position. Equa- 
tion (10) may be written 
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where x is the location of the center of gravity in 
fraction of chord as measured from the leading edge. 
Since vu s vf. 



Propellers usually operate near a Mach number of 
one, and the compressibility correction therefore becomes 
extremely important. As yet, there is no accurate knowl- 
edge concerning the compressibility correction for the 
flutter velocity near the velocity of sound. An approxi- 
mate compressibility correction for the subsonic range 
from reference ij. is 


M 


2 

c 




where M c is the Mach number corresponding to flutter 
speed in compressible flow and M is the Mach number 
corresponding to flutter speed in incompressible flow. 
The flutter speed corrected for compressibility is tenta 
tively calculated in the appendix and is indicated in 


figure 6 by the vertical line at 



The choice of the radius of the representative sec- 
tion is open to some question. Since the velocity varies 
approximately as the radius, this choice is rather impor- 
tant. It has been customary to use the section at three- 
fourths semispan as the representative section for wings. 
Because of the velocity distribution on propellers, the 
representative section was taken at the 0.8-radius station. 
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CONCLUDING REMARKS 


It has been shown that the stall flutter speed of a 
propeller Is in general very much lower than the calculated 
classical flutter speed. 

The classical flutter speed may be attained only if 
the propeller operates at the ideal angle of zero twist. 

The ideal angle of zero twist depends on the moment coef- 
ficient of the section, and the corresponding lift coeffi- 
cient has been given by a simple relation. 

It is desirable to have the design angle equal to 
the ideal angle of attack in order that the speed at which 
flutter occurs may be higher. The design angle should 
therefore be equal both to the ideal angle and to the 
Ideal angle of zero twist. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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. . APPENDIX 

SAMPLE CALCULATION OF FLUTTER SPEED CORRECTED 
FOR COMPRESSIBILITY FOR PROPELLER A 


Propeller section characteristics at 0.8-radius stations 
Type - Clark Y flat-bottom 12-percent- thick airfoil 
x = O.ljlj- 
r a 2 = 0.2i|. 

Specific gravity = 0.5 
1 

K ‘ k5 
~ = 0.053 


Propeller characteristics; 
R = 1.87 ft 
w a = 2n ( 535 ) 



f( 0.2lj.)(ij.5)(0.25) 
0.14}. - 0.25 


= 772 fps 


= ( 0.092) ( 2ir) (555 ) 
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Correction for compressibility: 



v f = ( 0 . 6 12 ) ( 1120 ) = 685 f*Ps 
c 



This value represents the flutter speed corrected 
for compressibility. (See fig. 6.) 
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Figure 1.- Force* on a representative propeller section. 
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Figure. 4«- Data for propeller A. 




Tip bending flutter amplitude, in 




Figure 5.- Bending and torsion amplitudes and flutter frequency of 
propeller B as function of flutter speed. Torsional frequency in 
still air, 340 cycles per second. 
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Fig. 6 
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Figure 6.- Flutter speed q/q cr as function of 
Ql u for propeller A. qAlcr 3 0*79 is the 
calculated flutter speed corrected for 
compressibility. 
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Fig. 7 
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Figure 7.- Twist of propeller tip as function of 
for propeller A. Tip speed constant at q/q cr 


- 0.37. 
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